Comparing GW+DMFT and LDA+DMFT for the testbed material SrVQ 3 



(N 

o 

(N 

> 
O 



o3 



i 

o 
o 



> 



X 



C. Taranto 1 , M. Kaltak 2 , N. Parragh 1 - 3 , G. Sangiovanni 1 ' 3 , G. Kresse 2 , A. Toschi 1 , and K. Held 1 
1 Institute for Solid State Physics, Vienna University of Technology, 1040 Vienna, Austria 
2 University of Vienna, Faculty of Physics and Center for Computational 
Materials Science, Sensengasse 8/12, A-1090 Vienna, Austria and 

3 Institut fur Theoretische Physik und Astrophysik, 
Universitat Wiirzburg, Am Hubland, D-97074 Wiirzburg, Germany 
(Dated: Version 0, November 7, 2012) 

We have implemented the GW+dynamical mean field theory (DMFT) approach in the Vienna ab 
initio simulation package. Employing the interaction values obtained from the locally unscreened 
random phase approximation (RPA), we compare GVF+DMFT and LDA+DMFT against each 
other and against experiment for SrVC>3. We observed a partial compensation of stronger electronic 
correlations due to the reduced GW bandwidth and weaker correlations due to a larger screening of 
the RPA interaction, so that the obtained spectra are quite similar and well agree with experiment. 
Noteworthily, the GW+DMFT better reproduces the position of the lower Hubbard side band. 
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I. INTRODUCTION 

The local density approximation (LDA) plus dynam- 
ical mean field theory (DMFT) approach-^— has been 
a significant step forward for calculating materials with 
strong electronic correlations. This is, because — on top 
of the LDA — -DMFT^i includes a major part of the 
electronic correlations: the local ones. In recent years, 
LDA+DMFT has been applied successfully to many ma- 
terials and correlated electron phenomena, ranging from 
transition metals and their oxides to rare earth and their 
alloys, for reviews see Refs. 0IE 

For truly parameter-free ab initio calculations, how- 
ever, two severe shortcomings persist: (i) the screened 
Coulomb interaction is usually considered to be an ad- 
justable parameter in LDA+DMFT and (ii) the so-called 
double counting problem (i.e., it is difficult to deter- 
mine the electronic correlations already accounted for at 
the LDA level). These shortcomings are intimately con- 
nected with the fact that the density functional nature 
of LDA does not match, with the many-body, Feynman- 
diagram structures of DMFT. Hence, it is not clear what 
correlations are already included on the LDA level or how 
to express these as a self-energy to avoid a double count- 
ing with DMFT correlations. These problems can be 
mitigated, but not solved, by constrained LDA (cLDA) 
calculations^—, which can be exploited to extract two 
independent parameters: interaction and double count- 
ing correctionj2iii 

A conceptionally preferable and better defined many- 
body approach is achieved if one substitutes LDA by the 
so-called GW approximatio n 12 ' 13 . Since its proposition 
by Biermann et al^, the development and application 
of such a GVF+DMFT scheme for actual applications 
has been tedious. This is reflected in the number of 
LDA+DMFT calculations for actual materials, which is 
of the order of a few hundred, compared to a single 15 
GVF+DMFT calculation for Niii, despite many advan- 
tages of GW +DMFT, such as the possible rigorous defi- 



nition in terms of Feynman diagrams and the avoidance 
of introducing ad hoc parameters for the Coulomb in- 
teraction and double counting corrections. The reason 
for this imbalance is twofold. First, since the GW ap- 
proach is computationally fairly demanding and com- 
plex, mature GW programs were missing in the past. 
Second, the GW+DMFT scheme is considerably more 
involved than LDA+DMFT, in particular, if calculations 
are done self-consistcntly and with a frequency dependent 
(screened) Coulomb interaction. Indeed, these concepts 
are presently tested on the model leve l 16 ' 20 . Let us also 
note in this context that a frequency dependent interac- 
tion has been employed on top of an LDA bandstructure 
for BaFe 2 As2^ and SrVO^ 



In this paper, we present results of our GW+DMFT 
implementation in the Vienna ab initio simulation pack- 
age (VASP) for SrV0 3 and compare the GW+DMFT 
results to those of LDA+DMFT— as well as photoemis- 
sion spectroscopy^ 2 -. We find the GW+DMFT spectra to 
be quite similar to that of LDA+DMFT due to a partial 
cancellation of two effects: the reduced GW bandwidth 
in comparison to LDA and the weaker screened Coulomb 
interaction. An important difference, however, is the po- 
sition of the lower Hubbard band which in GW+DMFT 
better agrees with experiment. To mimic the frequency 
dependence of the Coulomb interaction, which we have 
not included, we also performed GW+DMFT calcula- 
tions for a Bose factor ansataiS reduced bandwidth. The 
obtained spectra are rather different from GW+DMFT 
without Bose factor and LDA+DMFT. 



The paper is organized as follows: In Section UH we 
discuss the method and implementation. In Section IlIIl 
we compare LDA+DMFT and GW+DMFT self energies 
and spectral functions. A comparison to photoemission 
experiments is provided in Section IIVI and a conclusion 
in Section fVl 
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II. METHOD 

Let us briefly outline the relevant methodological as- 
pects. Starting point of our calculation is the GW im- 
plementation within (VASP). 30 Specifically, we first per- 
formed Kohn Sham density functional theory calcula- 
tions using the local density approximation for SrV03 
at the LDA lattice constant of a = 3.78 A and deter- 
mined the Kohn Sham one-electron orbitals 4>nk and one- 
electron energies e„k- The position of the GW quasipar- 
ticle peaks were calculated by solving the linear equation 

E nk = e "k + Z nU x 

Re[(<f> n k\T +V n - e + V H + £(e„ k )|<£ 

where T is the one-electron kinetic energy operator and 
V n - e and Vjj are the nuclear-electron potential and the 
Hartree-potential, respectively. £ is the GqWq self en- 
ergy, and Z„k is the renormalization factor evaluated 
at the Kohn-Sham eigenvalues ! 29 ' 30 The original Kohn 
Sham orbitals are maintained at this step. The Kohn 
Sham orbitals expressed in the projector augmented wave 
(PAW) basis are then projected onto maximally local- 
ized Wannier functions^ 3 - using the Wannier90 code 2^ 
To construct an effective low-energy Hamiltonian for the 
tig vanadium orbitals, we follow Faleev, van Schilfgaarde 
and Kotani and approximate the frequency dependent 
GqWq self-energy by an Hermitian operator H that re- 
produces the position of the quasiparticle peaks of the 
original self-energy exactly; 42 ' 43 

H mnM = i[(^, k |E*(^) + E(^ k P )Kk>]. (2) 

In practice, for the present calculations, we have applied 
the slightly more involved procedure to derive an Her- 
mitian approximation outlined in Ref. l3lL although this 
yields essentially an identical Hermitian operator H mn ^. 
Furthermore the off diagonal components are found to be 
negligible small, and henceforth disregarded. The final 
Hermitian and k-point dependent operator H is trans- 
formed to the Wannier basis and passed on to the DMFT 
code, where it is used to construct the k-dependent self 
energy by adding the local DMFT self energy. 

Fig. [1] shows the obtained GqWq bandstructure, which 
for the t 2g vanadium target bands is about 0.7 eV nar- 
rower than for the LDA. The oxygen p band (below 
— 2eV) is shifted downwards by 0.5 eV compared to the 
LDA, whereas the vanadium e g bands (located about 
1.5 eV above the Fermi- level) are slightly shifted upwards 
by 0.2 eV. In the LDA, the top most vanadium ti g band 
at the M point is slightly above the lowest e g band at 
the r point, whereas the GqWq correction opens a gap 
between the ti g and e g states. 

For brevity we will refer to the subsequent DMFT cal- 
culations as GW+DMFT, although it should be kept in 
mind that no self-consistency is performed in the GW 
part and furthermore the GW self energy is approxi- 
mated by an Hermitian (frequency in-dependent) non- 
diagonal operator. Subtracting the local part of this 



Hemritian operator (to avoid a double counting) yields 
for the degenerate ti g orbitals a constant shift. This pro- 
cedure allows to maintain the structure and outline of the 
common DFT-DMFT scheme and can be easily adopted 
in any DMFT code. However it neglects lifetime broaden- 
ing introduced by other bands not treated on the DMFT 
level. 

Within this Wannier basis, we also calculate the 
screened Coulomb interaction using the random phase 
approximation (RPA). As described in Ref. |35j for an ac- 
curate estimate of the interaction value to be used in 
DMFT (C/ DMFT ), only the local screening processes of 
the ti g target bands of SrVC>3 are disregarded since only 
these are considered later on in DMFT. This approach 35 
is similar to the constrained RPA (cRPA) 35,37 , with 
the difference being that in cRPA also non-local screen- 
ing processes of the ti g target bands are disregarded 
which are not included in DMFT. We carefully compare 
GW+DMFT with LDA+DMFT calculations and exper- 
iment. In both cases we use (frequency-independent) 
interactions obtained from this locally unscreened RPA 
and cLDA. The Kanamori interaction parameters as 
derived from the locally unscreened RPA are: intra- 
orbital Coulomb repulsion {7 DMFT = 3.44 eV; inter- 
orbital Coulomb repulsion £7 DMFT = 2.49 eV; Hund's ex- 
change and pair hopping amplitude J DMFT = 0.46 eV 4^ 
These values are, for SrV03, almost identical to the 
cRPA. 35 In cLDA, on the other hand, somewhat larger 
interaction parameters were obtained and are employed 
by us for the corresponding calculations: £/ cLDA = 
5.05 eV, U cLDA = 3.55 eV, J cLDA = 0.75 eV^ 2 . 

For the subsequent DMFT calculation, we employ 
the Wurzburg-Wien w2dynamics code 3 - 9 -, based on the 
hybridization-expansion variant^ - of the continuous-time 
quantum Monte Carlo method (CT-QMC) 4 !. This al- 
gorithm is particularly fast since it employs additional 
quantum numbers for a rotationally-invariant Kanamori 
interaction 3 - 9 -. The maximum entropy method is em- 
ployed for the analytical continuation of the imaginary 
time and (Matsubara) frequency CT-QMC data to real 
frequencies 4 ^. 

All our calculations are without self-consistency, which 
is to some extend justified for SrV03. Since the three 
tig bands of SrVC>3 are degenerate, DMFT does not 
change the charge density of the low-energy ti g mani- 
fold and hence self-consistency effects are expected to be 
small for LDA+DMFT. This is, in principle, different 
for GW+DMFT. Here, the frequency dependence of the 
DMFT self energy might yield some feedback already for 
a simplified Faleev, van Schilfgaarde and Kotani quasi- 
particle self-consistenc y 42 ' 43 . Finally, we also test the 
Bose factor ansatai 9 - which renormalizes the GW band- 
width by a renormalization factor Zb = 0.7— for mim- 
icking the frequency dependence of the (partially un- 
screened) RPA interaction. 
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FIG. 1: (Color online) Upper panel: GqWo quasiparticle 
bands (red, grey) in comparison to LDA (black). The Fermi 
level sets our zero of energy and is marked as a line. Lower 
panel: Wannier projected ti B bandstructure from GoWq (red, 
grey) and LDA (black). The t2 a target bands bandwidth is 
reduced by ~ 0.7 eV in GW. 



FIG. 2: (Color online) Comparison of the imaginary part of 
the DMFT self energies £ vs. (Matsubara) frequency cj for 
SrV03 at inverse temperature = 40 eV -1 as computed in 
five different ways: employing GW and LDA Wannier bands, 



the locally unscreened RPA interaction U 



= 2.49 eV 



and the cLDA U = 3.55 e V, as well as the Bose factor 
renormalization of Zb = 0.7.— 



III. RESULTS 

For analyzing the differences between GVF+DMFT 
and LDA+DMFT we analyze and compare in the fol- 
lowing five different calculations: 

1. LDA+DMFT@t7 cLDA (conventional LDA+DMFT 
calculation with the cLDA interaction (7 cLDA = 
3.55 eV). 

2. LDA+DMFT@(7 DMFT (LDA+DMFT calculation 
but with the locally unscreened RPA interaction 
t/ DMFT = 2.49eV). 

3. Giy+_DMFT@(7 DMFT (GVF+DMFT calculation 
with L> DMFT = 2.49eV) 

4. GPF+DMFT@(7 cLDA (GVF+DMFT calculation 
but with (7 cLDA = 3.55 eV) 

5. GPF+DMFT@(7 DMFT ,Z B = 0.7 (as 3. but with a 
Bose renormalization factor Zb) 

Let us first turn to the imaginary part of the local self 
energy which is shown as a function of (Matsubara) fre- 
quency in Fig. [2] The self energy yields a first im- 
pression how strong the electronic correlations are in 
the various calculations. The LDA+DMFT@t7 DMFT self 
energy is the least correlated one, somewhat less corre- 
lated than LDA+DMFT@t7 cLDA due to the smaller lo- 
cally unscreened Coulomb interaction (JJ BMFT = 2.49 eV 



< 3.55 eV = C/ cLDA ). For the same reason also the 
GVF+DMFT@(7 DMFT self energy is less correlated than 
that of a GW+DMFT@U cLDA calculation. 

If we compare LDA+DMFT and GW+BMFT on the 
other hand, the LDA+DMFT self energy is less corre- 
lated than the GVF+DMFT one, if the Coulomb interac- 
tion is kept the same. This is due to the 0.7 eV smaller 
GW ^-bandwidth in comparison to LDA. This observa- 
tion also reflects in the DMFT quasiparticle renormaliza- 
tion factors Z which were obtained from a fourth-order 
fit to the lowest four Matsubara frequencies, see Table 
HI Note that there is an additional GW renormalization 
factor reducing the bandwidth in comparison to LDA. 

However, the effect of the smaller GW bandwidth par- 
tially compensates with the smaller £7 DMFT interaction 
strength. Altogether this yields rather similar self ener- 
gies of the standard_ approaches: LDA+DMFT@(7 cLDA 
and GT4 7 +DMFT@t7 DMFT , see lower panel of Fig. [2 
This also reflects in very similar renormalization factors 
in Table H Z = 0.51 vs. Z = 0.57, which both agree well 
with experimental estimates of 0.5-0.6 j 22 i 45 

Since one important difference is the strength of the 
interaction, it is worthwhile recalling that jj r ' MFT is de- 
fined as the local interaction strength at low frequen- 
cies. While this value is almost constant within the 
range of the ^-bandwidth, it is much larger at larger 
energies, exceeding 10 eV. It has been recently argued 
and shown in model calculations^ that the stronger 
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LDA+DMFT@£/ cLDA 


0.51 


0.004 


0.013 


0.009 


LDA+DMFT@£7 DMFT 


0.67 


0.007 


0.016 


0.013 


GW / +DMFT@[7 DMFT 


0.57 


0.005 


0.014 


0.010 


GW+DMFT@i7 cLDA 


0.39 


0.003 


0.010 


0.007 


GW / +DMFT@(7 DMFT ,-Zs = 0.7 


0.36 


0.003 


0.009 


0.006 


experiment 22 ' 45 


~0.5-0.6 









TABLE I: (Color online) DMFT quasiparticle renormaliza- 
tion factors Z from the five different calculations at in- 
verse temperature /3 = 40 eV -1 . Also shown are the 
pairwise double occupations within the same orbital dmtra 
and between different orbitals with the same rioter an d op- 
posite spin 4J. W . The "standard" LDA+DMFT@£/ cLDA 
and GW+DMFT@[/ DMFT calculations are similarly corre- 
lated and well agree with experiment. Using the cLDA 
interaction([7 cLDA ) for GW+DMFT or the locally unscreened 
RPA ((/ DMFT ) for LDA+DMFT yields a too strongly and 
too weakly correlated solution in comparison to experiment, 
respectively. Note that GW+DMFT becomes even more 
strongly correlated, if the Bose renormalization factor is in- 
cluded. 

frequency-dependence of the screened Coulomb interac- 
tion at high energies is of relevance and can be mim- 
icked by a Bose factor renormalization of the GW band- 
width. The latter has been determined as Zb = 0.7 for 
SrVC>3. We have tried to take this into account in the 
Giy+DMFT@C7 DMFT ,Z B = 0.7 calculation. Due to the 
additional bandwidth renormalization, this calculation is 
very different from all others and yields the largest quasi- 
particle renormalization, i.e., Z = 0.36 is smallest. 

Next, we compare the k- integrated spectrum in Fig. 
[3J At low-frequency we find the same trends as for the 
self-energy results: The "standard" GW+DMFT and 
LDA+DMFT at J7 DMFT and U cLDA , respectively, yield 
a rather similar spectrum. In particular the quasiparticle 
peak has a similar weight and shape. However, a differ- 
ence that we will come back to later on is found at larger 
frequencies: The GW+DMFT Hubbard bands are closer 
to the Fermi level in comparison to LDA+DMFT. If we 
perform GW+DMFT and L_DA+DMFT_at the "wrong" 
interaction strength (i.e., £/ cLDA and {7 DMFT ; respec- 
tively), we obtain a noticeable stronger and weaker cor- 
related solution. This trend is also reflected in the double 
occupations presented in Table HI Finally, as in the case 
of the self-energy, the GW+DMFT@^ DMFT ,Z B = 0.7 so- 
lution is much more strongly correlated, with Hubbard 
side bands at much lower energies. 



IV. COMPARISON TO PHOTOEMISSION 
SPECTROSCOPY 

An obvious question is whether LDA+DMFT or 
GW+DMFT yields "better" results. This question is dif- 
ficult to answer and for the time being we resort to a com- 
parison with experimental photoemission spectroscopy 



(PES) 22 . However, one should be well aware of the lim- 
itations of such a comparison. On the theory side, the 
involved approximations common to the calculations, as 
e.g. neglecting non-local correlations beyond the DMFT 
and GW level, or further effects, such as the electron- 
phonon coupling or the photoemission matrix elements, 
might bias the theoretical result in one way or the other. 
On the experimental side, care is in place, as well. The 
PES results considerably improved in the last years due 
to better photon sources. Furthermore, in Ref. [22| an 
oxygen p-background has been subtracted, which by con- 
struction removes all spectral weight below the region 
identified as the lower Hubbard band. 

Fig. [4] compares the proposed LDA+DMFT and 
GW+DMFT (with and without Bose renormalization) 
with PES experiment. To this end, the theoretical re- 
sults have been multiplied with the Fermi function at 
the experimental temperature of 20K and broadened by 
the experimental resolution of 0.1 eV. The height of the 
PES spectrum has been fixed so that its integral yields 
1, i.e., accommodates 1 ^g-electron, as in theory. 

The GW+DMFT@[7 DMFT and LDA+DMFT@l7 cLDA 
have a quite similar quasiparticle peak, which also well 
agrees with experiment, as it was already indicated by the 
quasiparticle renormalization factor. A noteworthy dif- 
ference is the position of the lower Hubbard band which 

is at -2eV for LDA+DMFT@C/ cLDA and 1.6 eV for 

GW+DMFT® [/ DMFT . The latter is in agreement with 
experiment and a result of the reduced GW band width. 
Let us note that the sharpness and height of the lower 
Hubbard band very much depends on the maximum en- 
tropy method, which tends to overestimate the broaden- 
ing of the high-energy spectral features. Hence, only the 
position and weight is a reliable result of the calculation. 

As we have already seen, the Bose-factor renor- 
malized GVF+DMFT@C7 DMFT ,Z B = 0.7_ calculation 
is distinct from both, GW+DMFT@[/ DMFT and 
LDA+DMFT@C/ cLDA . It is also different from experi- 
ment with a much more narrow quasiparticle peak and a 
lower Hubbard band much closer to the Fermi level. A 
similar difference between static U on the one side and 
frequency dependent U was reported in Ref. 0. A dif- 
ference of this magnitude is hence to be expected. In the 
course of finalizing this paper, we became aware of Ref. 
0, in which Tomczak et al. report a GW+DMFT cal- 
culation with the full frequency dependence of the cRPA 
interaction for SrV03 obtaining good agreement with ex- 
periment as well. 



V. CONCLUSION 

We have carried out a careful comparison of 
LDA+DMFT, G W +DMFT and experiment for the 
case of SrV03, which is often considered to be a "bench- 
mark" material for new methods. To this end, the 
LDA or Go Wo quasiparticle bandstructure was projected 
onto maximally localized Wannier orbitals for the ti g 
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FIG. 3: (Color online) Spectral function for SrVC>3 computed in five different ways as in Fig. [2] At lower temperatures the 
central peak gets only slightly sharper and higher, although the temperature effects from /3 = 25 to 40 eV -1 are small. 
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FIG. 4: (Color online) Comparison of LDA+DMFT@t/ cLDA , 
GVK+DMFT@t/ DMFT (without and with Bose renormaliza- 
tion Zb = 0.7) and experiment. The position of the lower 
Hubbard band is better reproduced in GVK+DMFT whereas 
the central peak is similar in LDA+DMFT and GW+DMFT. 
The Bose renormalization GW+T>MFT differs considerably 
(photoemission spectra reproduced from Ref. |22T). 



bands. For these in turn correlation effects have been 
calculated on the DMFT level. If we take the lo- 
cally unscreened RPA interaction (or the similar cRPA 
one) for the GW+DMFT and the cLDA interaction 
for LDA+DMFT, the two approaches yield rather sim- 
ilar self energies and spectral functions at the Fermi 
level. These also agree rather well with photoemission 
spectroscopy. A noteworthy difference between these 
two calculation is found, however, for the position of 
the lower Hubbard band, which is better reproduced in 
GVF+DMFT. Similar spectra were also obtained by Tom- 
czak et al^ using a GIF+DMFT calculation including 
the frequency dependence of the interaction. 

We thank R. Arita, S. Bicrmann, K. Nakamura, Y. No- 
mura, M. Imada, S. Sakai, and J. Tomzcak for many hel- 
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Austrian Science Fund (FWF) through project ID 1597 
(as part of the DFG Research Unit FOR 1346) and the 
SFB ViCom F4103. Calculations have been done on the 
Vienna Scientific Cluster (VSC). 
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